Proper Dirac Quantization of Free Particle on D-Dimensional Sphere 
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We show that an unambiguous and correct quantization 
of the second-class constrained system of a free particle on a 
sphere in D dimensions is possible only by converting the con- 
straints to abelian gauge constraints, which are of first class in 
Dirac's classification scheme. The energy spectrum is equal to 
that of a pure Laplace-Beltrami operator with no additional 
constant arising from the curvature of the sphere. A quan- 
tization of Dirac's modified Poisson brackets for second-class 
constraints is also possible and unique, but must be rejected 
since the resulting energy spectrum is physically incorrect. 

1. Quantization of a free point particle in curved space 
is a long-standing and controversial problem in quantum 
mechanics. Dirac has emphasized that canonical quanti- 
zation rules are consistent only in a Cartesian reference 
frame ||]. Attempts to generalize these rules to curved 
space run into the notorious operator- ordering problem 
of momentum and coordinates, making the Hamilto- 
nian operator non-unique. Podolsky avoided this 
problem by postulating that the Laplacian in the free 
Schrodinger operator H = —h A/2 should be replaced by 
the Laplace-Beltrami operator Alb = g~^^^d^g^^'^g^^du, 
where 9^ = d/dq'^ are the partial derivatives with respect 
to the Z?-dimensional curved-space coordinates, and g 
is the determinant of the metric 5^1/ (g). This postu- 
late has generally been accepted as being correct since it 
yields, for a _D-dimensional sphere of radius R embedded 
in a Z? + 1 -dimensional Cartesian space with coordinates 
2;% an energy L^/2R^. Here La = —ipi{La)ijXj with 
Pi = —ihd/dxi are the unique quantum-mechanical dif- 
ferential operator representation of the D{D + l)/2 gen- 
erators La of the rotation group SO{D + 1) in flat space. 
If we take a to label the index pairs ij with i < j, then 
La are the matrices iL^j)ki = i{5ik5ji - SikSji), whose 
nonzero commutation rules are 



[Lij,Lik] = -iLjk (no sum over i) 



(1) 



The canonical commutation rules [pi,Xj] — —ihSij trans- 
fer the Lie algebra (|l|) to the operators Lij/h. The square 
of the total angular momentum L^ is the Casimir opera- 
tor of the orthogonal group S0(D+1), with the eigenval- 
ues Z(Z + D - 1), ? = 0, 1, 2, . . . 0. Thus, by quantizing 
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classical angular momentum rather than canonical vari- 
ables (i.e., by putting hats on L's rather than on p's) 
operator-ordering problems are avoided, making this en- 
ergy spectrum the most plausible one pi . 

Doubts on the correctness of this spectrum have been 
raised by DeWitt in his first attempt to quantize the 
system by a straightforward generalization of Feynman's 
time-sliced path integral to curvilinear coordinates. He 
found an extra energy proportional to the Riemannian 
scalar curvature R of in the Schrodinger operator: 



H=- — Ai^B+an'R 



(2) 
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with a proportionality constant a — 1/24. Various suc- 
cessors have presented modifications of DeWitt's proce- 
dure leading to other proportionality factors a = 1/12 
§ and a = 1/8 §. 

For the above sphere, the extra term produces an extra 
constant energy equal to ah D{D — 1)/ R^ which would 
contradict the previous results. At fist sight, this con- 
tradiction seems to be rather irrelevant. Common exper- 
iments are only capable of detecting energy differences ^ 
in which this constant drops out. Cosmology, however, is 
sensitive to an additive constant, which would change the 
gravitational energy of interstellar rotating two-atomic 
gases. Apart from such somewhat esoteric applications, 
a correct quantization of this simplest non-Euclidean sys- 
tem is certainly of fundamental theoretical interest. 

It is therefore important to find physical systems for 
which the ^-term is not a constant, so that it can dis- 
tort the energy spectrum of the Schrodinger operator 
(^). Such a system has recently been found, although 
in a somewhat indirect way. When solving the path in- 
tegral of the hydrogen atom, a two-step nonholonomic 
mapping is needed which passes through an intermediate 
non-Euclidean space ^ ,|l^ . The existence of an extra ^- 
tcrni in (0) would modify significantly the known level 
spacings in the hydrogen spectrum. Thus experiment 
eliminates an extra curvature scalar, in agreement with 
Podolsky and with the quantization of angular momen- 
tum. The successful procedure was generalized into a 
simple nonholonomic mapping principle, which allows to 
map classical and quantum-mechanical laws in fiat space 
into correct laws in curved space |9|-fil[. 

The absence of an ^-term in H) is therefore a cru- 
cial test for any quantum theory in curved spaces. It is 
the purpose of this note to show that this absence can 
be established for a sphere also within the operator ap- 
proach to quantum mechanics, after suitably preparing 
the description for an application of Dirac's theory of 
constrained systems. 



This application is, however, not straight-forward. A 
ZJ-dimensional sphere is most simply described by em- 
bedding it into a, D + I -dimensional Cartesian x-space 
via a constraint x^ — E? — 0. Within Dirac's classifi- 
cation scheme H, such a constraint is of second class. 
We shall demonstrate that Dirac's quantization rules for 
such systems produce wrong energy levels. A correct 
quantization becomes possible only by making use of a 
recently developed conversion [^ of second-class con- 
straints to first-class constraints, which instead of con- 
figuration space restrict the quantum-mechanical Hilbert 
space B. This conversion requires an extension of the 
phase space of the initial Cartesian system to a larger 
auxiliary Cartesian phase space, where the correct quan- 
tization is known. The operators representing the first- 
class constraints are generators of gauge transformations, 
and the physical states are all found by going into the 
gauge invariant subspace of the Hilbert space. 

Constraints associated with gauge symmetry have first 
been mastered in quantum electrodynamics (Coulomb's 
law) , and are now a standard tool in the quantization of 
gauge theories Q . 

2. We begin by quantizing the system with the help of 
Dirac's theory for second-class constraints. A free point 
particle with a Hamiltonian H = p^/2 moving in a flat 
Euclidean D + 1-dimensional is restricted to the surface 
of a Z?-dimensional sphere by the primary constraint in 
configuration space [E[ 



ip^=x^ - R^ ^0 



(3) 



The dynamical consistency condition ipi = {H,Lpi} = 
leads to an additional secondary constraint in phase space 



IP2 = (X,p) = 



(4) 



It expresses the fact that a motion on a sphere has 
no radial component. Although the canonical variables 
are Cartesian, the canonical quantization is not applica- 
ble, since the constraints (H), (P cannot be enforced for 
the associated operators — the conditions (pi = ip2 = 
would be in conflict with the commutation relation 
[(^1,(^2] = ^ihx-^ = 2ihR^ ^ 0. To resolve this diffi- 
culty, Dirac replaced the Poisson symplectic structure 
{xi,pj} — Sij by the so-called Dirac brackets 



{A, B}d = {A, B} - {A, ipa}[A-'r{ipb, B} 



(5) 



where Aab is a matrix A^h = {ipa,fb} formed from all 
primary and secondary constraints ipa = 0. This matrix 
is assumed to be non-degenerate, a defining property of 
second-class constraints pi ) . The Dirac brackets provide 
us with an antisymmetric operation which is distributive 
and associative, i.e., which satisfies Leibnitz rule and Ja- 
cobi identity, thus forming a symplectic structure which 
is as good as Poisson's. 

The removal of the inconsistency is ensured by the au- 
tomatic property 



{A, ipa}D = , for any A and ipa 



(6) 



Hamiltonian equations of motion generated by the 
Dirac bracket A = {A, H}d coincide with those gen- 
erated by the original Poisson bracket { , } on the sur- 
face of constraints (pa = 0. Thus, Dirac brackets pro- 
duce classically the same equations of motion as Pois- 
son brackets. Quantization proceeds by the replacement 
{, Id— >— i[,]/7i. The property m) allows us to replace 
the constraints by operator equations ipa — without the 
earlier contradictions. Moreover, since the constraint op- 
erators (fa commute with any other operator, they can 
be given any c-number value, for instance zero. For the 
surface of the sphere, the quantized Dirac brackets (||) 
read 



yXi , Xj\ 







[xj,pk] = ih ( Sjk T2- 



[PjiPk] = i^—{PjXk -PkXj) 



(7) 
(8) 

(9) 



The operator ordering problem occurring in the right- 
hand side of (y) is uniquely resolved under the condition 
that the algebra (R)-(ro satisfies the Jacobi identity. The 
operator x^ commutes with all canonical operators and 
is therefore a c-number, which can be set equal R^. 

To find the spectrum of the Hamiltonian H = |5^/2, 
we make use of the identity 



ylja^ )i\^aX )j X^Xj 



Jij 



(10) 



Inserting this into p — piSijPj, we obtain the Hamilto- 
nian 



H=^2Ll + \p\l 



(11) 



where La = —ip.i{La)ijXj as before and p,. = {x,p)/R 
is a radial momentum operator. Note that the operator 
(O) is determined without operator-ordering ambigui- 
ties. The identity (|o|) can equally well be inserted to 
the left and the right of the momentum operators piPj , 
with a unique result. 

The radial operator pr commutes with all canonical 
operators and is therefore some complex number c. It has 
necessarily an imaginary part, due to the obvious relation 
Pr—pl = ihD / R implied by (ra). Thus we may decompose 
Pr ^ c ^ ihD/2R + Cr, where Cr is an arbitrary real 
number. The constant Cr is determined by expressing 
Pr in terms of the operator of the second constraint (H). 
The ordering of p and x in it is undetermined, so that 
ip2 niay be written as {pr +Pr)/'2 + i?i'y with an arbitrary 
imaginary part 7. Inserting the above constant for pr 
and setting ip2 equal to zero is solved uniquely by 7 = 
and c,- = 0. 

It is now easy to find the spectrum of the first term in 
the Hamiltonian operator (O). The modified canonical 
rules (|3)^(§) transfer the Lie algebra (|l|) to the operators 
La — —ipLaX in just the same way as ordinary canonical 



commutation rules, so that the energy levels of angular 
momentum I are 



E^^^m + D-D + ^^n^D^ 



(12) 



The additional constant energy has a different in- 
dependence than the ^-term proposed by DeWitt and 
others Jt],^. In particular, it is nonzero even for a parti- 
cle on a circle {D = 2), where nobody ever expected such 
a term. 

Thus the modified symplectic structure proposed by 
Dirac, although esthetically appealing and yielding a 
unique result for a particle on the surface of a sphere, 
must be rejected on physical grounds as contradicting 
the established quantization via angular momentum op- 
erators and nonholonomic mapping principle. 

3. In second-class constrained systems, not every 
phase space variable can be made an operator. Dirac's 
new symplectic structure (||) represented by (0)-(^ ac- 
counts for this fact via its degeneracy in the embedding 
phase space. This permits the canonical variables of fluc- 
tuations transverse to the manifold on which the particle 
moves to remain c-numbers. There is, however, a defect 
in Dirac's procedure: Although physical excitations of 
the transverse degree of freedom are eliminated, the sys- 
tem maintains a memory of the forbidden motion by a 
nonzero c- number valued transverse energy p^pr/ 2. From 
the physical point of view, the existence of such a mem- 
ory must be rejected. After all, the embedding space is 
only an artifact for the introduction Euclidean canonical 
commutation rules. It does not belong to the manifold 
where the particle moves. An alternative approach must 
therefore be found where a memory of the embedding 
space is absent. 

Such an alternative is offered by gauge theories. In 
these, equations of motion exhibit symmetries which are 
local in time with the consequence that the dynamics 
of some degrees of freedom is not specified by the equa- 
tions of motion. Such local symmetries are gauge sym- 
metries, and the undetermined degrees of freedom corre- 
spond to pure gauge configurations. Upon quantization, 
non-physical degrees of freedom are removed by restrict- 
ing the Hilbert space to a physical subspace formed by 
gauge invariant states. 

Thus, rather than restricting the particle motion to the 
surface of a sphere via constraints as above, we consider a 
free motion in a Cartesian coordinate system, and impose 
the condition that the physical states in Hilbert space 
are invariant under arbitrary time-dependent rescalings 
of the radial size of the system. The transverse momen- 
tum pr can be made a generator of gauge transforma- 
tions, and we may require p}.Pr to be zero for physical 
states. 

To achieve this goal, we invoke the method of abelian 
conversion that allows one to transform a second-class 
constrained system into an abelian gauge theory p3[ . 
Dynamics of physical (gauge invariant) degrees of free- 
dom in the effective abelian theory is the same as dy- 



namics of physical degrees of freedom in the original 
second-class constrained system. In general, the abelian 
conversion proceeds as follows. Given a set of second- 
class constraints ipa {a = 1,2, ... , 2M), one extends the 
original phase space by extra independent canonical vari- 
ables Qa,Pa (a = 1,2,...,M). The extended phase 
space is equipped with the canonical symplectic structure 
{Qq, Pji} = Sap and {xi,pj} = Sij (all other brackets are 
zero) . With the help of abelian conversion, quantum dy- 
namics on a manifold can be formulated independently 
of the parametrization of the manifold [E3 . 

An equivalent set of abelian first-class constraints CTq = 
is constructed in such a way that it satisfies 



{aa,<Jb} ^ 0. 



(13) 



This amounts to solving first-order differential equations 
with the boundary condition aa{x,p,P — 0,Q — 0) — 
(pa{x,p). Given the new constraints aa = 0, the original 
system Hamiltonian H{x,p) is converted into a Hamilto- 
nian on the extended phase space H(x,p, P, Q) by solving 
the equation 



{il,aa} = 



(14) 



with the boundary condition H{x,p,Q — 0,P = 0) ~ 
H{x,p). The extrema of the associated extended action 
S"^^^ = J dt{piXi+PaQa—H — XaCra) determine cquatious 
of motion in the extended phase space. These depend on 
2M arbitrary functions of time Xa(t) as a manifestation 
of the gauge freedom. There exists a choice for Xa{t) such 
that the auxiliary phase space variables Pa (t) and Qa (t) 
vanish at all times, whereas pi{t) and Xi{t) solve the orig- 
inal equations of motion of the second-class constrained 
system |13{ . 

Applying this procedure to our particular second-class 
constraints (0) and (0) yields 



CTl 



(pi+ P , a2 ^ (p2 + 2x Q , 



while the extended Hamiltonian assumes the form 



H 



1 



+ 



Ll 



2 \ai+R^ cri+i?2 



(15) 



(16) 



where La = —ipLaX are the classical components of the 
angular momentum (for a = ij. La = XiPj — XjPi). The 
extended phase space variables p, x, P, Q are Cartesian 
and satisfy the standard Poisson bracket relations. They 
can be turned into hermitian operators in the usual way 
by the replacement { , } — > — i[ , ]/h 

Physical states are invariant with respect to transfor- 
mations generated by ai^2'- x — > e~^^x and Q ^ Q — £,i, 
where ^1^2 are parameters of the gauge transformations 
uM . The first is an arbitrary time-dependent rescaling of 
the size of the system, which is geometrically equivalent 
to the initial restriction of the motion to the surface of the 
sphere, while the second implies that the auxiliary degree 



of freedom is a pure gauge. The first-class constraints re- 
strict the physical Hilbert space to the gauge-invariant 
sector by the Dirac conditions M 



CTl,2^phys = 



(17) 



The general solution has the form \&phys = f{Q, x'^)'^{^): 
where f{x,Q) is some fixed function, whereas ^(ri) are 
wave functions on the /^-sphere. In the physical Hilbert 
space, we can set cri_2 to zero in the Hamilton operator 
(O). Thus we find the energy values 



Ei = ^il + D-1), 



(18) 



rather than dT3) . There is no additional constant energy, 
in agreement with P,B,P|- 

This result is obtained without ordering problems. Al- 
though ordering ambiguities are not absent altogether, 
they do not affect the final result since they occur only 
in the quantization of the gauge generator (T2, where they 
modify only the explicit form of the physically irrelevant 
function f{x'^, Q), but not the physical Hilbert space de- 
scribed by the wave function ^(il), nor the spectrum 
(p^). In fact, in the simple system at hand the ordering 
ambiguities in tT2 produce only a multiplicative renormal- 
ization of the physical states *phys- 

4. The above quantization of a particle on a sphere 
via abelian conversion is of course applicable to arbitrary 
homogeneous spaces. But what about arbitrary mani- 
folds? When attempting a straight-forward generaliza- 
tion we run once more into operator-ordering problems 
for the extended Hamiltonian (|l6|), and these require a 
new strategy. The quantum Hamiltonian H = p^/2 has 
no ordering ambiguity and is again adopted as a starting 
point. Let n{x) be a unit vector normal to the man- 
ifold in the embedding space (if this space has a di- 
mension higher than D + 1, more normal vectors are 
needed to specify transverse directions). The condition 
Pn^phys = in{x) , p)'^ phys = rcmovcs the transverse 
motion and offers itself as a generator of gauge transfor- 
mations. This, however, is not consistent for two reasons. 
First, since n(x) depends on position, [iJ,p„] is nonzero, 
so that the free-particle Hamiltonian is not gauge invari- 
ant. Second, the operator pn is not hermitian making 
finite gauge transformations non-unitary. 

The first problem can be resolved via an abelian con- 
version method performed immediately at the quantum 
level flJ. Here one uses operator versions of second- 
class constraints tpi = F{x) and Lp2 — Pn to restrict 
the motion to a manifold specified by F{x) = 0, the 
forbidden direction being specified by the normal vector 
n{x) = dF{x)/\dF{x)\. Then one extends the system by 
two extra canonical operators Q and P obeying standard 
Heisenberg commutation relation, and commuting with x 
and p. The conversion of the second-class constraint op- 
erators ipa is enforced by solving equations (03) and (llj) 
with Poisson brackets replaced by commutators. The 



abelian gau ge g enerators have the same basic structure 
as those in (15), only that the i-dependent factor of Q 
in (72 depends now on F{x). This construction solves the 
first problem of finding a gauge invariant Hamiltonian 
operator. 

The hermiticity problem for the generator (72 must be 
solved in a way that dynamics of physical degrees of free- 
dom does not depend on the embedding procedure. For 
this we observe that the operator cr2^2 is just as good as 
(72 itself to eliminate excitations of the transverse modes. 
It has the advantage of being hermitian. In addition, it 
commutes with (7i, so the abelian gauge algebra (^ is 
retained by such a choice of the second generator. Thus 
we modify the conversion method by taking plj)n as the 
constraint operator Lp2 rather than p„. At the classical 
level, the new constraint p^ = is certainly equivalent 
to Pn = 0. But at the operator level, it is superior by 
generating unitary gauge transformations. With the new 
generator, the conversion equation ( |l4[) gives rise to a 
new Hamiltonian operator in the extended Cartesian co- 
ordinate system. 

It is not hard to verify that the physical states ( p^ ) 
have the form e^^^/^"^{x) where pl^Pn^ — 0, that is, the 
kinetic energy of the transverse motion is strictly zero for 
physical states with our choice of the gauge generators. 
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